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Abstract. The embedded ensembles were introduced by Mon and French as 
physically more plausible stochastic models of many-body systems governed by one- 
and two-body interactions than provided by standard random-matrix theory. We 
review several approaches aimed at determining the spectral density, the spectral 
fluctuation properties, and the ergodic properties of these ensembles: moments 
methods, numerical simulations, the replica trick, the eigenvector decomposition of the 
matrix of second moments and supersymmetry, the binary correlation approximation, 
and the study of correlations between matrix elements. 
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1. Introduction 

Canonical random-matrix theory (RMT) as introduced by Wigner |3Sj considers 
ensembles of random matrices classified [THj by their symmetries. In the asymptotic 
limit of infinite matrix dimension, N — > oo, this theory makes a number of remarkable 
predictions. The average spectrum (the spectral density) has the shape of a semicircle, 
the spectral fluctuations are universal (i.e., are under very weak conditions independent 
of the weight function used to define the ensemble, are parameter-free, and dependent 
only on the symmetries of the ensemble), and the results are ergodic (so that for most 
members of the ensemble the spectral average of any observable coincides with the 
ensemble average of that same observable). For a recent review of the field, we refer to 
reference |22j . 

Nevertheless, early applications of canonical RMT to nuclear spectra encountered 
objections or, at least, questions. This is because in RMT, the number of independent 
random variables (i.e., of uncorrelated matrix elements) is of order N 2 . Put differently, 
in RMT every state in Hilbert space is connected to every other such state by a matrix 
element which does not vanish and is assumed to be an independent random variable. 
But with the success of the nuclear shell model, it had become clear that nuclei are 
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governed by one- and two-body interactions (the mean field and the residual two-body 
interaction). In a representation of Hilbert space spanned by Slater determinants of 
single-particle wave functions (solutions of the mean-field equations), the residual two- 
body interaction has non-vanishing matrix elements only between determinants which 
differ in the occupation numbers of not more than two single-particle states, and most 
matrix elements of the residual two-body interaction, therefore, vanish. Even if one 
is prepared to accept a stochastic approach and considers the non-vanishing two-body 
matrix elements of the residual interaction as uncorrelated random variables, the ensuing 
stochastic model is very different from canonical RMT: The number of independent 
random variables is generically small compared to the dimension of a typical shell- 
model matrix. Why — so the question — can canonical RMT then serve as a model for 
the description of fluctuation properties of nuclear spectra? 

French and Wong [TBI CHj and Bohigas and Flores [HI CD] approached this question 
with the help of numerical simulations. These authors replaced the actual two-body 
matrix elements of the nuclear shell model by random variables (this, in essence, 
defines the two-body random ensemble (TBRE)) and studied the resulting spectral 
fluctuation properties. They concluded that the fluctuation properties were consistent 
with predictions of canonical RMT. 

Because of the complexities of angular-momentum and spin coupling and the 
exclusion principle, the TBRE is not amenable to analytical investigations. Thus, the 
introduction of the fc-body embedded ensembles (EGE(fc)) by Mon and French [32*] 
might be considered an essential step towards an analytical understanding of the 
numerical results obtained by French and Wong [TS1 [JJ"] and Bohigas and Flores j*2 EI] 
and, thereby, of the stochastic properties of nuclear spectra. Since then, the question 
has gained much wider significance: Canonical RMT has successfully been used [22] i n 
such diverse many-body systems as atoms, molecules, and quantum dots. All of these 
systems share with nuclei the feature that they are known to be governed essentially by 
effective one- and two-body forces. Thus, the embedded ensembles may be viewed as 
the generic models for stochasticity in many-body systems. 

The embedded ensembles dispose of all the complexities of the couplings of angular 
momentum and spin while retaining the symmetries imposed by the exclusion principle. 
One considers m Fermions in I degenerate single-particle states interacting via a random 
/c-body interaction. The single-particle states carry no further quantum numbers like 
spin or angular momentum. In contrast to shell-model calculations (where the single- 
particle states are usually not degenerate) degeneracy is here assumed in order to focus 
attention on the results of a stochastic A;-body interaction. The interaction is not 
restricted to k — 2 (although this remains the physically most interesting case). This 
is done in order to obtain an understanding of the transition from the case k = 2 to 
the case k = m which, as it turns out, is equivalent to canonical RMT. Just as in 
RMT, the embedded ensembles can be bestowed with unitary, orthogonal or symplectic 
symmetry, and one distinguishes accordingly EGXJE(k), EGOE(/c), and EGSE(/c). As 
in canonical RMT, one is interested in universal properties. Therefore, one considers 
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the limit of infinite matrix dimension. This limit is attained by letting the number 
of single-particle states go to infinity, I — > oo. EGE(/c) is, thus, a generic model for 
stochasticity in many-body systems governed by A;-body interactions. Needless to say, 
EGE(fc) can and has been generalized to the case of Bosons. 

Central questions in the theory of the embedded ensembles then are: (i) What is 
the shape of the spectral density? (ii) What are the spectral fluctuation properties? 
(iii) Are these properties universal? (iv) Are the spectra ergodic? 

In this Review, we present the status of the field in the light of these questions. 
We start with a definition of EGE(fc). Then, we order the material by the analytical or 
numerical technique which has been used, rather than by the properties of EGE which 
have been addressed. 



2. Definitions: The /c-body Embedded Ensembles of Fermions and Bosons 

We define the /c-body embedded ensembles for Fermions and for Bosons for the cases 
of unitary and of orthogonal symmetry. We do not consider symplectic symmetry 
throughout this paper. 

We consider a set of / degenerate single-particle states labelled \ j) , with j = 1, . . . , I. 
The associated creation and annihilation operators are denoted by aj and a,j in the case 
of Fermions and fej and bj in the case of Bosons. These operators obey the usual 
(anti)commutation relations. Then, \j) = aj|0) or \j) = bj\0), with |0) the vacuum 
state. 

To define the /c-body interaction, it is useful to introduce operators which create 
a normalized state with k < I Fermions or k Bosons from the vacuum. In the case of 
Fermions, these are written as 

k 

rf,a = ^Lw* = n4> (!) 

s=l 

where ji < j 2 < ■ ■ ■ < jk- The label k defines the number of single-particle creation 
operators, while a stands for the set of j's. The corresponding annihilation operators 
are given by 

lk,a = (4 J ■ (2) 

In the case of Bosons, the A;-body operators have the form 

k 

Xk,a = x] 1 j 2 ,...j k = N a l[bl, (3) 

s=l 

where ji < j 2 < . . . < jk- Here, N a is a normalization factor which guarantees that 
the state xl JO) has norm unity. For every set of equal indices j containing n elements 
occurring in a, N a contains the factor 1/y/nJ.. 

The random A;-body interaction for Fermions is given by 

V k = J2 V k;a^l a ^k, 7 - (4) 

«7 
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The coefficients ffc ;a7 are random variables with a Gaussian probability distribution, 
mean value zero, and a common second moment. The form an ensemble of random 
operators. In the case of unitary symmetry, the Vk m are complex. The only non- 
vanishing second moment is 

The bar denotes the ensemble average. The Kronecker delta's stand for the string 
Sjij'djzj' ■ ■ ■■ In the case of orthogonal symmetry, the ffc ;a7 are real and obey 

Vk;a'yVk; a >y> = V 2 [5 aQ '5 7 y + <5 Q7 '5 7Q '] . (6) 

Following the usage in canonical RMT, we distinguish both cases by writing Vk(/3) 
with j3 — 1, 2 for orthogonal and unitary symmetry, respectively. The fc-body random 
interaction for Bosons is defined in complete analogy to equation (J3J. The parameter 
v 2 sets the scale of the energy. As long as we consider only a single EGE(fc), we can put 
v 2 = 1 without loss of generality. 

We consider a Hilbert space obtained by distributing m < I Fermions or m Bosons 
over the single-particle states \j). A complete set of basis states for Fermions (Bosons) 
is given by V'mJO) (by XmJ0)> respectively). The dimensions of these two spaces are 

N — ( ] for Fermions , 
\m J 

N B =( l + m ~ 1 ) for Bosons . (7) 
\ m J 

Here and in the sequel, we use the index B for Bosons while we typically work without 
such an index for the more frequent case of Fermions. We take m > k and consider 
the matrix of the fc-body interaction in that m-particle Hilbert space. The matrix 
elements (/i|\4|^) have the form 

W*\") = (olVwWJUo) , (8) 

and analogously for Bosons. Using equation (|3J) for 14, it is easily seen that for m — k, 
the matrices of the fc-body interaction coincide with the GUE or the GOE, depending 
on the symmetry chosen, both for Fermions and for Bosons. This is not the case for 
m > k. In this case, we speak of the fc-body embedded Gaussian ensembles (EGE(fc)) 
of random matrices in the m-particle Hilbert space. We use the notation EGUE(fc) and 
EGOE(/c) for the cases of unitary and orthogonal symmetry, respectively. Aside from 
the (unitary or orthogonal) symmetry, the embedded ensembles are defined in terms of 
the three parameters l,m,k. We recall that in canonical RMT, the limit iV — >• 00 of 
infinite matrix dimension yields universal results. For the EGE(fc), we proceed likewise 
by taking the limit I — > 00 which implies the limit of infinite matrix dimension. For 
Fermions, we have the obvious constraint that k < m < I. The limit / — > 00 with fixed 
k can be attained by keeping m fixed, or by keeping the ratio m/l < 1 fixed. Brody 
et al. define the dilute limit by taking / — > 00, m — > 00, m/l — > 0. For Bosons, 
m may be larger than I. Since EGE(fc) is constructed in such a way that EGE(m) 
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coincides with canonical RMT for both Fermions and Bosons, and since in the dilute 
limit the distinction between Fermions and Bosons becomes irrelevant, we expect that 
the spectral properties of EGE(A;) do not differ qualitatively for Fermions and Bosons 
except for the case where the number of m Bosons is close to or exceeds /, the dimension 
of the underlying single-particle space. This is indeed what is found. 

Canonical RMT is invariant under unitary or orthogonal transformations of Hilbert 
space. While this statement carries over trivially to the fc-body ensembles embedded 
into a Hilbert space spanned by k particles, that invariance is lost for m > k. The 
calculation of spectral properties of canonical RMT greatly benefits from such an 
invariance property. Thus, the determination of the spectral properties of the embedded 
ensembles is much more difficult than for canonical RMT. 

3. Moments Methods 

The first analytical approach to EGE(fe), developed by Mon and French [32], calculates 
moments of the matrix elements of and related quantities. It has led to insights 
into the shape of the EGE(fc) spectra (the spectral density), has led to the development 
of statistical spectroscopy in nuclei, and has given a limited answer to the question of 
ergodicity. These topics are reviewed in references fUI2E]- The approach has mainly 
been used for Fermions in the context of the nuclear shell model, although it can also be 
applied for Bosons as in the Interacting Boson Model for nuclei. We confine ourselves 
to Fermions. 

In a Hilbert space with m Fermions, the p th moment of is defined as M p (m) = 
(1/-/V) tr[(\4) p ]. It is obvious that all odd moments vanish. The even moments can 
be evaluated using Wick contraction of the Vfc's. The moment is then the sum over all 
patterns of pairwise contractions. Assigning the same letter to pairs of Wick-contracted 
V^'s, the fourth moment is, for instance, given by M 4 (m) = (1/N) tr {^lA^A^B^B^ + 
AkBkAkBk) where we have used the cyclic invariance of the trace. For I — > oo, the 
low moments can be worked out using combinatorial arguments [32] • Knowledge of the 
moments yields information on the average spectral density p(E). This is because of 
the relation 



where the Ei are the eigenvalues of Vk- 

Canonical RMT is obtained for k = m and predicts the spectral density to have 
semicircular shape. On the other hand, in the dilute limit, all Wick-contracted pairs 
of 14 become independent of each other, the moment M 2p (m) is given by M 2p (m) = 
(2p — l)!!(M 2 (m)) p , and the spectral density is Gaussian. This shows that for fixed 
k and for m monotonically increasing from m = k, the spectral density undergoes a 
transition from the semicircular to the Gaussian shape. For k — 1, the transition can be 
worked out easily because the level density for m > 1 is the convolution of m semicircles. 
Already for m = 9, the spectral shape is very close to Gaussian. This argument fails for 




(9) 
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k > 2. In reference jS2], the even moments M 2p (m) are worked out for p = 1, 2, 3, 4. The 
values obtained are consistent with the gradual transition from semicircular to Gaussian 
shape of the spectral density. 

In nuclear theory, sums over all final states of the strength of some transition 
operator O play an important role. Examples are single-particle transfer strengths or 
the Gamow-Teller transition strength. Such transition strength sums can be worked 
out, too, using the moments method. The reason is that K = O^O can be expressed 
in terms of creation and annihilation operators, and terms like tr[(Vfc) p .K'] can again be 
calculated using combinatorial arguments, especially in the dilute limit. This method 
has been very successfully applied in the framework of the nuclear shell model. In many 
cases, EGOE(/c) results obtained in this way show remarkable agreement with results 
of detailed shell-model calculations. We do not dwell on this important issue because 
an extensive review has appeared quite recently [28J. 

Using the moments method to calculate the spectral density or transition strength 
sums, we always use the ensemble average. Is it justified to compare such averages 
with properties of a given physical system (as we do when comparing with shell-model 
calculations)? While a given physical system obviously does not permit averaging over 
the (physically non-existent and fictitious) ensemble, it does permit the determination 
of spectral averages. For instance, by grouping neighbouring levels into bins, we can 
determine the mean spectral density. This poses the question: Are spectral averages 
(taken over a specific system) and ensemble averages (over the EGE(fc)) identical? This 
problem is referred to as the problem of ergodicity, in analogy to the well-known ergodic 
problem in classical statistical mechanics concerning the equality of phase-space average 
and the long-time average over a single trajectory. 

For ensemble averaging to be relevant for an individual system, the property under 
study must vary little from member to member of the ensemble. In the context of the 
moments method, the question of ergodicity is answered by calculating the variances 

Sg)(m) = (l/iV)tr(\4 P ) WmWT) ~ NNMVf) ■ (10) 

For canonical random-matrix theory and in a normalization where the individual 
moments are finite in the limit N — ► oo, the variances vanish asymptotically: 
The canonical ensembles are ergodic as far as the spectral density is concerned. For 
the embedded ensembles, the calculation of the variances in equation ([TO)) involves 
contractions between pairs of operators located in different traces as well as in the 
same trace. The treatment of the first kind of pairs requires combinatorial techniques 
beyond the ones developed for the moments [22]. Mon and French could show that for 
EGOE(/c), £^(m) vanishes asymptotically for I — > oo. Thus, they could demonstrate 
ergodicity for the spectral density. Their argument can easily be extended to EGUE(fc). 

Unfortunately, the moments method, while very useful for the spectral density and 
for statistical spectroscopy, is only of limited use for calculating spectral fluctuations. 
For canonical RMT, Mon and French consider density fluctuations as excitations of the 
semicircle and use this approach to study spectral GOE fluctuations. They succeed in 
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re-deriving the logarithmic dependence of the spectral stiffness on the length of the 
energy interval. However, the method apparently does not permit the study of short- 
scale fluctuations like the nearest-neighbour spacing distribution. Mon and French apply 
the same method to EGOE(A;) with m ^> k. They show that long-range fluctuations of 
the spectrum are large. As a consequence, neither the position of a given (numbered) 
eigenvalue nor the position of the centre of the spectrum or the variance are ergodic. We 
return to this point in Section The short-range fluctuations, however, which yield 
information about spectral statistics, pose difficulties. To quote from reference [EE] : 
"...there is no real theory yet for EGOE fluctuations, the gap being one of the most 
significant ones in the entire subject." 

4. Numerical Results 

Insight into the spectral fluctuation properties of the embedded ensembles may be gained 
from numerical simulations. Such simulations have been performed from the early days 
of the field. Here we review separately the simulations for Fermionic and for Bosonic 
systems. 

4-1. Fermionic Systems 

Prior to the introduction of EGE(fc) by Mon and French [32], numerical simulations 
used the two-body random ensemble (TBRE). This ensemble was introduced in the 
context of nuclear physics where standard shell-model techniques were available to 
calculate spectra of given spin and parity (and isospin if relevant) in terms of a fixed 
two-body interaction. The latter was replaced by an ensemble of two-body matrix 
elements with a Gaussian distribution, zero mean value, and a common second moment. 
The antisymmetrized (Fermionic) random two-body matrix elements thus define the 
residual interaction among m particles distributed over I degenerate single-particle 
states. Inasmuch as nuclei are governed by one- and two-body forces, the TBRE is 
clearly a more suitable stochastic model than canonical RMT. The complexities arising 
from the angular-momentum and spin couplings, from correlations between many-body 
matrix elements induced by the two-body character of the interaction, and from Pauli's 
exclusion principle made it impossible to treat the TBRE analytically. This led French 
and Wong fSJ EH] and Bohigas and Flores [HI EI] to resort to numerical simulations 
using matrices of rather small dimensions (< 61). The main question addressed in these 
simulations was: Are the spectral fluctuation properties of the TBRE the same as those 
of the GOE? 

The main results of these early works are: (i) The spectral density of the TBRE 
is Gaussian, in agreement with the results of shell-model calculations with realistic 
interactions, (ii) The spectral density displays a transition from Gaussian to semicircular 
shape as the rank of the interaction k is increased, (iii) Unfolding the spectrum of each 
member of the ensemble separately yields good agreement of the fluctuation properties of 
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the TBRE with those predicted by canonical RMT. (iv) The fluctuations of the position 
of the first eigenvalue of the TBRE are stronger than those for the GOE. (v) The TBRE 
is neither stationary nor ergodic. These results hold also for simplified versions of the 
TBRE, which do not involve, or at least do so only partially, the complexity of the 
angular momentum algebra. 

Early results on the p th nearest neighbour spacing distribution displayed systematic 
deviations between the TBRE and the GOE JT] while experimental data showed 
agreement with GOE predictions. This seemed to suggest the existence of an effective 
many-body interaction in nuclei. It was soon realized, however, that the deviations 
were related to the method of analysis of the spectra and caused by the non-ergodic 
properties of the TBRE [3*3 120] ■ More precisely, the deviations were due to the difference 
between ensemble averaging and spectral averaging. Ensemble averaging uses a single 
staircase function (the integrated level density) to unfold the spectra of all members 
of the ensemble. This function is obtained by superposing the levels of all members 
of the ensemble. Early numerical results [TBI E] (and later analytical studies |2HE2]) 
had shown that in the dilute limit (I ^> m ^> k, I —>■ oo, m — > oo with m/l — > 0), the 
mean level density is Gaussian. This spectral shape was accordingly used for ensemble 
unfolding. On the other hand, spectral averaging uses a separate staircase function for 
each member of the ensemble. It was found that the deviations from GOE behaviour 
obtained by ensemble averaging were related to the variance of the spacing distribution; 
an analytical estimate of this effect, linking ensemble and spectral average, was also 
obtained [20] E2]- A first systematic comparison of spectral and (a corrected) ensemble 
averaging, using matrices of rather small dimensions, was carried out in reference j!3j . 
This comparison showed remarkable agreement between the TBRE results and the GOE 
predictions. The agreement invalidated the earlier conclusion that nuclear interactions 
are mainly of many-body type, and has left open the question of the rank of the effective 
nuclear interaction I14j. 

Unfolding the spectrum is technically a subtle task. It can considerably influence 
the results on spectral fluctuations. Attention was therefore devoted to a consistent way 
of constructing ensemble averages which yield correct spectral fluctuations f3J[2H]- This 
problem has recently again attracted attention ^3 EH!- m the dilute limit, the spectral 
density is Gaussian. However, knowledge of this fact is not sufficient for a proper 
unfolding of the spectra of the EGE(k) and of the TBRE as obtained in numerical 
simulations with finite-dimensional matrices. This is because for such matrices, the 
first and second moments of the spectrum are not ergodic (cf. Section f6.3|) . Hence, 
these moments must be normalized for every member of the ensemble. The sensitivity 
of the spectral fluctuations to these corrections was addressed in reference [TTj, where 
either the first or the second moment, or both, are corrected for every matrix of the 
ensemble. The results show that, in the center of the spectrum, the nearest-neighbour 
spacing distribution is not sensitive to such corrections. For the two-point function, 
ensemble unfolding with the Gaussian mean level density yields results that deviate 
clearly from predictions of canonical RMT. A slight correction results from adjusting 
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the width of the distribution for every member of the ensemble. A drastic improvement is 
obtained when the distribution for every member of the ensemble is recentered. Applying 
both corrections together yields results which are hardly distinguishable from those of 
canonical RMT. Once the low moments of the spectral distributions are adjusted so as 
to coincide for all members of the ensemble, a polynomial fit for the ensemble staircase 
function is usually implemented. This procedure yields a function with which the spectra 
can be meaningfully unfolded. 

A method for spectral unfolding which is somewhat more complicated but has a firm 
theoretical basis |3*2~1 ITi] for the embedded ensembles, uses a Gram-Charlier expansion. 
Starting from the fact that in the dilute limit the spectral density is Gaussian, one 
writes the spectral density p s (x) of an ensemble of matrices of finite dimension in the 
form [111 E3] 



Here, x — (E — E c )/a is the recentered energy in units of its standard deviation, pc 
is the Gaussian distribution and H n (x) is the Hermite polynomial of degree n. The 
Hermite polynomials describe the long-wavelength oscillations of p s (x) (long in terms 
of the mean level spacing) about the asymptotic value pg(x)- The coefficients c n are 
adjusted so as to minimize Arms, the overall root-mean-square error of the level-to- 
level deviations of the staircase function calculated from equation (JTTj) from the data. 
Explicit expressions for C3 and C4 are given in reference • One uses equation (fTTj) with 
some maximum value hq in the sum over n. This value should be sufficient to eliminate 
any secular trend in Arms while preserving the fluctuations of the spectrum. The value 
of n is usually determined by calculating Arms as a function of no, and by choosing 
the smallest value beyond which no significant corrections to Arms are found. For the 
EGE(/c) this method yields good agreement with the spectral fluctuation properties 
predicted by canonical RMT P3 ESI EH] • 

The spectral fluctuation properties at the edge of the spectrum have also received 
increased attention. These are interesting for a comparison of TBRE and EGOE(fc) 
results with the positions of low-lying nuclear states. Bohigas and Flores compared 
the properties of the low-lying part of the spectrum of the TBRE and of the GOE. 
They showed that the widths of the positions of individual eigenvalues were much 
larger for the TBRE than for the GOE. Cota et al. [T^j fitted the nearest-neighbour 
spacing distribution to a Brody distribution, P(s, u) = a(d; + l)s w exp[— as w+1 ] with a = 
+ 2) /(u + 1)]} W+ , obtaining for the Brody parameter the value id w 0.80 ± 0.05. 
More recent results by Flores et al show that the semi-Poisson distribution, 
P(s) = 4sexp[— 2s], gives a better fit than the Brody distribution, if the levels of 
each member of the ensemble are normalized first according to an individual Gaussian 
distribution (spectral unfolding). 

Numerical simulations using matrices of dimension ~ 3000 have repeatedly found 
agreement between the fluctuation properties of the EGE(fc) and the TBRE at the center 



Ps(x) = p G (x) 1 + 



n>3 



CnH n (x) 
n! 



(11) 
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of the spectrum on the one hand, and the predictions of canonical RMT on the other. 
Therefore, it is often taken for granted that a random k = 2 part in the Hamiltonian is 
sufficient to induce level fluctuations of canonical RMT type. This point of view is at 
odds with the fact that many-body systems with random few-body interactions (both 
for Fermions and Bosons) may possess a high degree of order in the low-lying part of 
the spectrum |2U 123 El ■ For instance, for the TBRE with all the angular momentum 
couplings taken into account explicitly, it was found that there is a statistical preference 
for J w = + ground states despite the fact that these states account only for a small 
fraction of Hilbert space [21]- Moreover, these J n = + ground states are separated by 
a gap from the lowest excited state. Other properties like localization in Fock space [2E] 
and odd-even binding effects [33J show that pairing effects are robust properties of two- 
body interactions. Such properties are usually understood in terms of Hamiltonians 
which involve some collective behaviour. All these results, based on Hamiltonians with 
random interactions, contradict the basic philosophy of and the predictions based on 
canonical RMT and suggest that the embedded ensembles may not always yield results 
which coincide with canonical RMT. 

4-2. Bosonic Systems 

Early analytical results by Kota and Potbhare |2Z| indicated a Gaussian spectral density 
both in the dilute and in the dense limits. The dense limit which exists only for Bosons, 
is defined as m — > oo for fixed I and k. Manfredi [SD], using matrices of dimension 364 
{I — 4, m — 11), compared for the first time the spectral fluctuation properties at the 
center of the spectrum of EGE(fc) for Bosons with the results of canonical RMT. Using 
spectral unfolding he concluded that there is no significant difference between the two 
ensembles. Patel et al [33] considered this problem in the dense limit, studying the case 
I = 5 and m = 10. These authors constructed the ensemble-averaged staircase function 
using a sixth order Gram-Charlier expansion (jllj) . They found excellent agreement with 
canonical RMT results for the nearest-neighbour spacing distribution (with a Brody 
parameter u = 0.95) and for the A 3 -statistic. These results led them to conclude that 
the embedded ensembles possess generically (for Fermions and Bosons) the spectral 
fluctuation properties of canonical RMT. 

Recently it was found analytically, however, that the embedded ensembles for 
Bosons are not ergodic in the dense limit [IJ Ej. More precisely, the fluctuations of 
the centroids and widths of the spectrum (in units of the average width) do not vanish 
in this limit, but attain constant values (cf. Subsection l(i.3j) . This fact implies that 
unfolding the spectra by ensemble averaging or by spectral averaging will yield different 
results. The Hilbert-space dimension for Bosons is given by equation (JJJ). It is clear 
that the effect can be best displayed numerically for 1 = 2. 

Figures ^ and El illustrate the non-ergodic behaviour for m = 3000 Bosons and 
the case of a two-body interaction (with 1512 members forming the ensemble) for 
the nearest-neighbour spacing distribution and for the A 3 -statistic, both taken in the 
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Figure 1. Nearest-neighbour spacing distribution P(s) obtained by (a) ensemble 
unfolding and (b) spectral unfolding. The dotted curves correspond to the Wigner 
surmise and the dashed ones to the Poisson distribution. We note the different vertical 
scales used in the frames. Taken from reference [2]. 




Figure 2. A3~statistics (solid lines) measured at the center of the spectrum after 
(a) ensemble unfolding and (b) spectral unfolding. The dotted curve corresponds to 
the GOE results, the dashed one to a Poissonian spectrum and the dotted-dashed one 
to a picket-fence spectrum. Taken from reference 



centre of the spectrum. Spectral unfolding (ensemble unfolding) was done by fitting 
a polynomial of degree 11 to the staircase function of each member of the ensemble 
(to the ensemble-averaged staircase function, respectively). For ensemble unfolding, 
the individual spectra were not recentered or rescaled since there is now no theoretical 
support for this procedure. In fact, rescaling and recentering the spectra would yield a 
non-Gaussian average level density. Interestingly, in the case of spectral unfolding the 
nearest-neighbour spacing distribution is dominated by a large peak centred at s = 1. 
This suggests that individual spectra have an almost constant level spacing, i.e., are 
close to spectra of picket-fence type. This is further supported by the large plateau 
observed in the plot for the A3-statistic. 

For the same case (I = 2, m = 3000) the structure of the eigenfunctions was 
also investigated [TJ E|. It was found that individual spectra consist generically of a 
superposition of independent sequences of levels with constant spacings as if there was 
some symmetry present in the system. When a new sequence appears, the old and the 
new sequence are interlaced, and the level spacing is no longer constant for a while. In 
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Figure 3. Intensities |ci„| 2 of a sequence of eigenvectors in the many-body basis 
in the vicinity of a kink in the staircase function. The two overlapping segments of 
nearly equidistant levels are easily distinguished by the structure of the corresponding 
eigenfunctions. Taken from reference PJ. 



the staircase function, this is seen as a sudden change in the average slope, 1.6. , cLS £L kink. 
The first few eigenstates of the new sequence are localized in Hilbert space, even though 
they may correspond to rather highly excited states of the spectrum. This is illustrated 
in figure El where the squares of the coefficients q„ of the eigenvector expansion in an 
ordered many-body basis \fi n ) are plotted (cf. reference |2j). 

These numerical results show that in the dense limit, the Bosonic EGE(fc) differ 
significantly from the canonical ensembles of RMT. In the case I — 2 it has been 
established that the difference is due to the exact integrability of the system in the 
semiclassical limit [Sj. Essentially, this case corresponds to the quantization of a classical 
two-degrees of freedom system which possesses two independent integrals of motion, 
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the energy and the number of particles. For arbitrary values of I, the corresponding 
classical system has I degrees of freedom. As I increases, chaotic trajectories appear in 
phase space, typically for energies around the center of the spectrum, and there is a 
corresponding transition in the spectral fluctuations of the quantum system. 

In concluding this Section we mention that the spectral properties at the edge of 
the spectrum have also been investigated for Bosons. Using a random interaction in 
the framework of the Interacting Boson Model, Bijker and Frank have demonstrated 
the emergence of collective motion in the low-lying part of the spectrum, and have 
studied the dependence of collective motion on the rank k of the interaction |K|. 
When the number of interacting Bosons m is large enough with respect to k, they find a 
preponderance of J n = + ground states as well as definite evidence for the appearance 
of vibrational and rotational structures. For m ~ k they do not obtain an indication of 
vibrational or rotational band structures. 

5. The Method of Replicated Variables 

Verbaarschot and Zirnbauer [37J used the replica trick developed in statistical mechanics 
for the study of spin glasses and of Anderson localization to investigate spectral 
properties of the GOE and of EGOE(A;). Using the replica trick, one writes the 
observable under study as the logarithmic derivative of a suitably chosen generating 
function Z, replaces logZ by lim n ^ [^ n — 1], evaluates Z n for positive integer values 
of n and takes the limit n —>■ of the resulting expression. Similar to what happens 
in the supersymmetry approach described in Section El below, averaging Z n yields in 
the exponent of Z n a quartic term in the integration variables which contains the 
matrix A^J. pa ((3) defined in equation |T2|) . The authors introduce matrix elements of 
the operators ip k which connect the m-particle states with the (m — A;)-particle states 
and perform the Hubbard-Stratonovich transformation in this mixed representation. 

They use the saddle-point approximation for the one-point function and find 
a semicircle for the spectral density. The moments method having shown that the 
semicircle does not apply for all values of (m, k), this points to the need to study the 
loop expansion. It is shown that there are terms in the loop expansion which correct 
the semicircle. However, it is not possible to evaluate all such corrections. 

The saddle-point approximation is then applied to the two-point function within 
the same mixed representation as used for the one-point function. The saddle-point 
solution allows for the existence of a Goldstone mode. This mode carries the variable 
s = (Ei — E 2 )/d where d is the mean level spacing. It is with respect to this variable 
that fluctuation properties of the spectrum are measured. The two-point function is 
evaluated using an expansion in inverse powers of s, the point being that the stiffness of 
the GOE spectrum relates to the occurrence of the factor s~ 2 in the average two-point 
function. Difficulties and, in fact, uncertainties arise from the need to integrate over 
the massive (i.e., the non-Goldstone) modes. With this proviso, it is found that the 
leading terms in the loop expansion combine to yield the value — (17/2)(l/(7r 2 s 2 )) for 
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the density-density correlation function. This result has to be compared with the GOE 
value — (l/(n 2 s 2 )). Led by early numerical results reviewed in Section 0J the authors 
believe that they should have obtained the GOE result and speculate that the difference 
is due to high-order terms in the loop expansion. They show that if the density-density 
correlation function is oc s~ 2 then stiffness of the spectrum is implied. 

Results obtained more recently [3 EJ [H |2| and reviewed in Sections |H1 and UJ cast 
some doubt on the assertion that the spectral fluctuation properties of EGE(A;) coincide, 
for all values of k, with those of the canonical ensembles, at least in the limit I — > oo. 
This would imply that higher-order loop corrections in the approach of reference [37] 
could modify the factor (17/2) in a fc-dependent fashion. 



6. The Matrix of Second Moments 



A novel approach to determine spectral properties of EGE(fc) was developed in 
references |S1 IH d E] • We describe this approach here for the case of Fermions. The 
case of Bosons can be treated in complete analogy. Starting point is the observation 
that the matrix elements (/x|Vfc|z/) of the stochastic operator Vk are themselves also 
random variables and have a Gaussian probability distribution with mean value zero. 
Therefore, all properties of the embedded ensembles are completely determined by the 
matrix A^J. p(T ((3) of second moments defined by 



lli , P M = {^\vMW)(p\yMW) ■ (12) 

The idea is to extract information on EGE(fc) by studying the properties of the matrix 
Affi ((3). As in the case of canonical random matrices, we use the labels (3 = 1 and 
(3 = 2 for the orthogonal and unitary ensemble, respectively. Performing the average 
yields 



WI^sI^KpI^I^W) + <W*#IV'a|p>] 



(13) 



In the last of equations (jlHj) we have used the reality of the matrix element (u\iplip a \p) . 



6.1. Duality. Eigenvector Expansion 

In the unitary case, there is a connection between the matrices A^{2) for the /c-body 
interaction and A^ m ~ kS) (2) for the (m — /c)-body interaction. The relation is referred to 
as duality and reads 

4tU( 2 ) = 4^ ) (2). (14) 

We note the difference in the sequence of indices {pupa} on the two sides of this equation. 
The proof of equation (|14j) rests on the fact that for every m-Fermion state \p) and for 
every operator ip k)0l with ipk,a\p) 7^ 0, there exists a uniquely defined operator ^(m-fc)^ 
such that i[>( m -k),'y' l (>k,a\p>) — |0), the vacuum state. We emphasize that the duality 
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relation has nothing to do with particle-hole symmetry. It applies likewise to Bosons 
and can be extended to the case (3 = 1. The duality relation obviously assigns a special 
role to the case 2k = m. This is the reason why spectral properties of EGE(&;) change 
at 2k = m. 

For 2m > I, the available single-particle states are more than half filled, and it is 
tempting to use particle-hole symmetry to simplify the algebra. Rewriting in terms 
of hole operators and bringing the hole creation operators up front, one creates a sum 
of fc'-body interactions with k' = 1, 2, . . . , k. This is why particle-hole symmetry is of 
limited use only, except for the dilute limit. 

A very useful tool in the analysis of spectral properties of EGE(fc) is the eigenvector 
expansion of the matrix A^ k \ We take the unitary case, (3 = 2, for Fermions and 
consider the matrix A^ k ' of second moments as a matrix in the product space {/if}. 
In this space, A^ k \2) is a Hermitean matrix and can, therefore, be diagonalized. The 
eigenvalue equation reads 

E 4tU 2 )^ a) = A (s) (fc)C^ . (15) 

pa 

The index s = 0, 1, . . . labels different eigenvalues and the index a labels degenerate 
eigenvectors. The eigenvectors are orthogonal. We choose the normalization 

T,C { ; u a) Cif=N6 st 6 ab . (16) 

Provided that the eigenvectors form a complete set, the matrix A^ k > can be expanded 
in the form 

-i m—k 

4tU( 2 ) = v E E A(s \k)C^C^ . (17) 

In writing equation (|17|) we have anticipated the fact that for s > (m—k), the eigenvalues 
A^ vanish. Knowledge of the expansion (j!7)l . i.e., of the eigenvectors and eigenvalues 
of A^ k \ makes it possible to calculate the low moments of V*.. Moreover, using the 
expansion (fTTj) in the framework of supersymmetry allows for the use of the Hubbard- 
Stratonovich transformation. The orthogonal case is treated analogously. Combining 
the results for (3 = 1, 2, we obtain 

-i m—k 

4SU(0) = E E A(s) W [c$>c£> + 8vC$>c£>\ . (is) 

ly s=0 a 

To appreciate the significance of the eigenvector expansion equation (JUJ), it is useful 
to recall the form of the second moment of the GUE Hamiltonian H^ u , 



A 



2 



Hp,aHpu = —jSfiufipo- ■ (19) 

We have put v 2 = 1 for EGE(fc) and, by analogy, put A = 1 for the GUE. The fact that 
GUE coincides with EGUE(m) then implies that for k = m, equation (|T7j) reduces to 
equation ()19j) . Put differently, the one non-vanishing eigenvalue of the GUE is A^°'(m), 
and the associated non-degenerate eigenvector is Cj$ = 5^ u . All other orthogonal 
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eigenvectors belong to eigenvalue zero. The comparison shows that the eigenvalue 
expansion equation (fTTj) is the natural generalization to EGUE(fc) of equation (fT^j) for 
the second moment of the GUE Hamiltonian. The Kronecker delta's in equation (j!9j) 
express the unitary invariance of GUE. The fact that for k < m the sum over s in 
equation (j!7)) extends up to (m — k), with eigenvectors Cjfif) which differ from Kronecker 
delta's, is due to the fact that Vk does not possess this unitary invariance for k < m. 
We remark in parenthesis that the eigenvalue in equation (fTTJj) obviously has value unity, 
while it turns out that for m = k, the eigenvalue A^(m) in equation (fTTj) has the value 
(J^j. The difference is due to the fact that when we use second quantization to write 

GUE = V m and work out A^(2), we have to count the number of holes which is (J^j. 
In the appendix of reference |SZ|, results closely related to equations (fT2*Jl to (fTHJ) were 
derived. 

It remains to determine the form of the eigenvectors and eigenvalues. The 
eigenvectors Cjf^' have the form 

c% a) = WlA>> , (20) 

where a stands for the set (a, 7). The form equation (}20|) applies for all s = 0, 1, . . . , m 
independently of the actual value of the eigenvalue A^ s \k). Whenever two or more 
single-particle indices in the sets a and 7 coincide, special attention is required 0] . For 
this reason, the number D^ s ' of orthonormal degenerate eigenvectors Cjf^ belonging to 
fixed s (and, thus, the dimension of the corresponding subspace of Hilbert space) are 
given by 

^•^Cf-f 1 X ■ (21) 




It follows that J2T=o = ( m ) , the dimension of the product space {fiv}, so that the 
eigenvectors form a complete set. The eigenvalues are given by 

AW(*)=f m - a V'- m + *- a ) . (22) 

We observe that A^(A;) = for s > (m — k). This fact limits the sums over s in 
equations ()17ll8j) . Using the eigenvector expansion in the duality relation, equation (JHJ), 
generates useful identities. 

For Bosons, the eigenvalue expansion equation (fTHj) applies likewise. The 
eigenvectors have the form given in equation fl2()D . with ip replaced by Xi an d are 
normalized as in equation (|16jh with iV replaced by iVg. The dimension D^) of the 
subspace spanned by degenerate eigenvectors characterized by s with s — 0, 1, . . . , m is 
given by = 1 and, for s > 1, by 



D$ = T ' : ' ) -("'"") • (23) 



I + s — 1 
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Again, we have Y^=qDb = Nb> showing that the eigenvectors form a complete set. 
The eigenvalues are given by 



The eigenvalues vanish for s > m — k, in keeping with the GUE result where for k = m 
the only non-vanishing eigenvalue belongs to s — 0. 

6.2. Group-Theoretical Aspects 

Duality and the eigenvector decomposition equation (JTHJl apply to both Fermions and 
Bosons. A group-theoretical analysis jSSj shows that both these results and much of 
the structure displayed in Subsection 16.11 apply much more widely and are rooted in 
symmetry properties of the embedded ensembles. The analysis also gives a deeper 
meaning to the concept "embedded ensemble". For simplicity, we confine ourselves to 
the unitary case although all results apply likewise to the orthogonal one. We denote by 
a] the creation operator for a Fermion or a Boson in single-particle state j and by ty ka 
the creation operator of a normalized state containing k particles (Fermions or Bosons). 
This normalized state need not be a Slater determinant or a product state (as is the case 
for the operators ipka and \ka defined in Section EJ) but may be a linear combination of 
such states. 

Three symmetry groups are relevant for the embedded ensembles. These are (i) the 
group SU(Z) of unitary transformations of the I degenerate single-particle states \ j) with 
j — 1, . . . , I; (ii) the group U(iVfc) of unitary transformations of the /c-body interaction 
where N k is the dimension of the Hilbert space containing k particles (Fermions or 
Bosons); (iii) the group U(iV m ) of unitary transformations of the Hilbert space containing 
m particles, with N m the dimension of this space. As we shall see, the group SU(/) 
governs the embedding, the group U(N k ) is obviously the symmetry group of EGUE(fc), 
and XJ(N m ) is the symmetry group of the GUE. 

An element u G SU(/) generates a unitary transformation of the single-particle 
states and a corresponding transformation T{u)a^T\u) of the creation operators aj. 
Under u, the operators ^>\ a transform according to the irreducible representation D£l(u), 



Here t is an integer which may take the values m, k, or (m — k). For Fermions (Bosons), 
the matrices D are in essence totally antisymmetrized (symmetrized, respectively) 
powers of u. 

In analogy with the familiar fractional-parentage technique, we can expand the m- 
particle states \ma) = ^FjnJO) into products of states containing k and (m — k) particles, 
respectively, 




(24) 



T{u)^ ta T\u) = Y.D^{u)y 



(25) 




(26) 
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The coefficients CJ"^ s are the coefficients of fractional parentage or, equivalently, 

J fc i J ( 771 — fa ) 

the Clebsch-Gordan coefficients for the coupling [(fkf(m-k))fm) °f t ne irreducible 
representations / fc and /( m _ fc ) to f m . 

The product ^jL^fca of operators which appears in the definition of the k- 
body embedded ensemble transforms according to the direct product of irreducible 
representations D^ k {u) and D^ k {u). The product can be reduced to a direct sum of 
irreducible representations D 9b (u). This defines a set of basic &;-particle interactions 
Bk(ba) which are Hermitean, transform according to the irreducible representation 
D 9b (u), and are given by vector-coupling \jjjL and ^j-s v i & the Clebsch-Gordan 

coefficient C 9 , h % r . Using the Wigner-Eckart theorem, we can write the matrix element 

fk-yf k ° 

(m , -f\B k (ba)\m5) as the product of a reduced matrix element (m\\B k (b)\\m) and the 
Clebsch-Gordan coefficient C 9 . ba -T r . 

The interaction V k can now be rewritten as a sum over the operators Bj.(ba). 
The coefficients are uncorrelated Gaussian distributed random variables. Using this 
form and the Wigner-Eckart theorem, one finds for the second moment of the matrix 
elements (m^\Vk\m8) an expression which is equivalent to the eigenvalue decomposition 
equation (|T5|). The role of the eigenvectors is played now by the Clebsch-Gordan 
coefficients, and the role of the eigenvalues is taken by the squares of the reduced 
matrix elements (m\\Bk(b)\\m) . The index b coincides with the summation index s 
in equation ()18|) . The duality relation also follows from the Wigner-Eckart theorem. 
Thus, all the relations derived explicitly for Fermions and Bosons are seen to follow from 
group-theoretical considerations. 

These insights allow for a generalization of the concept of an embedded ensemble 
which brings the underlying concepts to the fore most clearly. We consider an arbitrary 
compact simple Lie group G and two independent systems labelled k and (m — k) whose 
basic states \fkCt) and \f( m -k)l) transform according to the irreducible representations 
D fk (g) and D f ( m - h )(g) of the group G, with g e G. Now let us assume that a non-trivial 
interaction Vk of GUE type occurs only in system k. The embedding of this interaction 
into a space of different dimension is accomplished by projecting the product states 
\fk&) \f(m-k)l) onto the subspace of states which transform according to the irreducible 
representation D' m {g) which is contained in the direct product D^ k (g) x D^ m ~ k ){g). Let 
the associated projection operator be denoted by I(f m ). The embedded ensemble is 
then defined by the Hamiltonian I(f m )VkI(fm). This definition comprises the essence of 
the group-theoretical extension of the idea of an embedded ensemble. It is independent 
of the existence of Fermions and Bosons and relies only on group-theoretical concepts. 

Group-theoretical arguments can also be used to isolate that part of Vk which is 
invariant under U(N m ), and to investigate symmetries of the generating functional in 
the supersymmetry approach [S3]. The part of which is invariant under U(N m ) 
transforms either like the GUE or like a multiple of the unit operator. Both possibilities 
actually occur in the decomposition. The latter, taken by itself, would cause Poissonian 
level statistics. It carries particularly large weight for k -C m. 
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6.3. Moments ofV k 

The eigenvector expansion equation (|18j) and the orthogonality relations equation (|16|) 
make it possible to calculate low moments and variances of Vk without resorting to the 
dilute limit / 3> m, at least for (3 = 2. Three observables of interest are 



s ((l/N)trV k ((3)y 



(l/JV)tr(W)) s 



R = ((1/NM VM) 2 ) 2 " ( (l/N)tr(V k ((3)) 2 f 



(l/N)trV k (l3)) 



2 \2 



t= JWW, (27) 

( (1//V)trl4(,3)) 2 r 

The ratio S measures the fluctuations of the center of the spectrum in units of the 
average width of the spectrum. The ratio R measures the relative fluctuation of the 
width of the spectrum. The parameter k is the kurtosis, and Q marks the difference 
in spectral shape between the semicircle (Q = 0) and the Gaussian (Q = 1). Explicit 
values for S, R, Q are given in references [U 12], both for Fermions and for Bosons. The 
result for R was first obtained by French [20]. Both S and R vanish in the limit I — > oo. 
However, for fixed values for k and m/l, both ratios decrease very slowly (with inverse 
powers of the logarithm) with increasing dimension N of the matrices. This fact is 
at the root of the difficulty in obtaining reliable spectral information on EGE(A;) from 
numerical simulations. 

For the shape of the spectrum, the quantity Q is of primary interest. It is explicitly 
given 0] in terms of the eigenvalues A^ s '(k) and the dimensionalities D^ B > introduced 
above. In the limit I — > oo, Q vanishes if 2k > m with both k and m fixed, and with 
both k/m and m/l fixed, while Q attains a finite and non-zero value for 2k < m with 
both k and m fixed, and with both k and m/l fixed. This suggests that the transition 
from the semicircular to Gaussian shape takes place for 2k < m. (It would take the 
study of higher moments to make this conclusion unambiguous). As mentioned above, 
the critical role played by the value 2k = m is attributed to duality. 

For Bosons in the dense limit, the ratios 5* and R do not vanish as m — > oo. 
Thus, the fluctuations of the centroids and widths of individual spectra do not vanish 
asymptotically. Hence, the Bosonic ensembles are not ergodic in the dense limit m — > oo 
with k and / fixed. 



6.4- Super •symmetry Approach 

The eigenvector expansion equation (j!8)) makes it possible to apply the supersymmetry 
approach. To see this, we recall that when one uses supersymmetry to calculate 
properties of the GUE, equation (JTHJ) is extremely helpful. Indeed, after averaging 
the generating function over the ensemble, there appears in the exponent a term 
which is quartic in the commuting and anticommuting integration variables. Owing 
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to equation (|19|). this quartic term can be written as the square of a bilinear form, and 
the latter is handled with the help of the Hubbard-Stratonovich (HS) transformation. 
Because of equation (fT%|) . a similar situation arises for EGE(£;). The form of this equation 
implies that the quartic term in the integration variables obtained after ensemble- 
averaging, is a sum of squares of bilinear forms, and the HS transformation again can 
be used. There is a price, however: Whereas only a single graded matrix a is needed 
to carry out the HS transformation for the GUE, the EGE(fc) requires the introduction 
of as many cr-fields as there are independent eigenvectors to non-zero eigenvalues in 
equation ()18|) . This number rises steeply as m increases from k = m, see equation (}2*Tj) . 

The saddle-point approximation to the integration over the cr-fields yields for 
EGE(/c) a semicircle for the spectral density, and universal Wigner-Dyson spectral 
fluctuations, just like for the GUE. The loop expansion, generated by expanding 
the cr-fields around the saddle-point solution, serves as a test of the saddle-point 
approximation. For the GUE, all terms in the loop expansion vanish asymptotically 
for iV — > oo. For EGE(/c), it is not possible technically to go beyond the lowest non- 
vanishing term of the loop expansion. For the spectral density, one finds that for 2k > m, 
this term vanishes asymptotically for I — > oo. This fact reinforces the conclusion in 
subsection 16 . 31 that the spectral density has semicircular shape for 2k > m. For 2k < m, 
finite corrections appear the form of which is consistent with the value of Q found in 
subsection 16.31 For the spectral fluctuations (the two-point function), the situation 
is more ambiguous. The lowest-order loop correction vanishes asymptotically for all 
values of k albeit it does so ever more slowly as k decreases from k = m. This might 
be consistent with universal Wigner-Dyson level statistics for all values of k. We note, 
however, that the lowest-order loop correction also produces non-universal terms. Such 
terms might also arise in higher order and not vanish asymptotically. 

In summary, the supersymmetry approach yields useful information especially 
for 2k > m. It lends independent support to the conclusion that the change from 
semicircular to Gaussian spectral shape sets in as k decreases and passes through the 
value 2k = m. For the spectral fluctuations, it is consistent with but not necessarily in 
support of Wigner-Dyson statistics for all values of k. However, for the physically most 
interesting case k = 2, it yields very little information. 

7. Binary Correlation Approximation 

None of the analytical approaches described above has been able to yield definitive 
information on the spectral fluctuation properties of EGE(/c) for the physically 
interesting case k <C m. There exists, however, another approach first introduced by 
Mon and French for the calculation of the one-point function |32j . This is the binary 
correlation approximation which applies in the dilute limit fc <C m C I, m/l as 
m, I — ► oo [Tlj . The binary correlation approximation can be generalized to investigate 
the two-point function [3J |1] . We first review the calculation of the one-point function 
of the EGUE(/c) in the version of Verbaarschot and Zirnbauer |3*7| . and then present that 
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of the two-point function [HI E] . We note that in the dilute limit the distinction between 
Fermions and Bosons is physically irrelevant, both cases yielding the same results. For 
definiteness, we consider the case of Fermions. To be free of singularities, we choose in 
this Section v 2 = [A^(fc)] -1 which implies a bounded spectrum of unit width. 



7.1. One-Point Function 



Following Verbaarschot and Zirnbauer |HZj , we expand the one-point function in a power 
series in 14, 

1 



1 



1 



^ 00 

N ^ z 

iV p=0 z 



(2f 



where z = E — in, E is the energy and i] > is an infinitesimal increment. In the 
second equation in ()28|) . we have interchanged the summation and the average over 
the ensemble, and we have used the Gaussian distribution of Vk which implies that 
only even powers of Vk contribute. The ensemble average is carried out using Wick 
contractions. Each contracted pair of V^'s is evaluated in the dilute limit, i.e., replaced 
by v 2 hS°>(k) = 1. In the term of order p, there are (2p — 1)!! different ways of pairwise 
contracting the V^'s. The result is 



(2p 



1 00 



pi 



2z 2 



(2p)\ 



(29) 



p=0 " " p=0 - 

This expression is evaluated using the technique of Borel resummation. We use the 
identity n! = / °° e~ l t n dt for the factor (2p)\. We interchange the summation and the 
integration. The sum over p yields an exponential function, thus leading to 



dtexp 



2z 2 



Equation ()3(Jj) converges provided Re \z 2 ] < 0. Writing z 



(30) 

-i\z\e 1 ^, we see that this 



condition is fulfilled for — | < 



< |. The integral is evaluated as follows. First, we 
choose as contour the straight line that joins the origin and z. Putting t = re 1 ^, we 
rotate this contour so that it comes to lie on the real axis. A new (real) integration 
variable u = r/\z\ is introduced. The resulting integral converges for all z, and is 
evaluated in a straightforward manner, yielding j3Zj 

.,2-1 r £ z - 

(31) 



exp 



erfc 



V2. 



Taking the limit 77 — > and recalling that the mean level density is given by p(E — irj) 
7r _1 Im[(y((£' — irj)] we obtain the Gaussian form for the average level density 

2 



p(E) = (27r)- 1/2 exp 



(32) 
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7.2. Two-Point Function 



Our presentation differs slightly from that in reference [I] and takes account of recent 
developments |3E]- Again, we work in the dilute limit and consider the quantity 



«,(*,,)= »<*M*>1, (33) 
g{zi) ■ g{z 2 ) 

where zf = E ± +e/2 and zf = E ± —e/2, with E and e the mean value and the difference 
of the energy arguments of the two Green functions, respectively. An upper plus 
(minus) sign denotes an infinitesimal positive (negative) imaginary energy increment, 
respectively. We are interested in values of e which are of the order of the mean level 
spacing and, thus, small compared to unity (the width of the spectrum). Thus, we have 
approximately \z\\ ~ \z 2 \. In order to obtain the connected part of the density-density 
correlator, z\ and z 2 must have imaginary energy increments of opposite signs. 

Proceeding as in the case of the one-point function, we expand R 2 (z\, z 2 ) in powers 
of Vfe. Using equation (j2Hj) we have 

OO OO J 



g(z 1 )g(z 2 ) = — £ £ _ w tr[(^)i'] tr[(V k y) . (34) 

iV p =0q=0 z l z 2 

The Gaussian distribution of the V k implies that p + q must be even. Using Wick 
contraction we obtain two types of contributions. The members of the contracted pair 
either occur under the same trace, or involve both traces. We refer to the latter as 
cross-contractions. Reordering terms, we obtain 



i oo oo 1 7 \ — 

^i)^ 2 ) = T^E E 2p+s+ i 5^ tr[OT tr[(U fc ) 2 «+1 • (35) 

ly s=0 p,q=0 Z l Z 2 V 7 s 

In equation f(35j) . the index s counts the number of cross-contracted pairs. There are 
( 2P s" S ) (^s^) different ways of cross-contracting the V^'s, while there are (2p— 1)!! (2g — 
1)!! ways of pairwise contracting the remaining Vjt's. For the latter, we use the binary 
correlation approximation. Thus, each one of these contractions yields a factor unity, 
irrespective of the position where the contracted V^'s appear in each of the traces. The 
result is 

7 w x 1 ^ ^ (2p-l)!!(2g+l)!! ^2p + ^ /2g + g\ 

iV s=0p,?=0 Z l Z 2 V 6 / V 6 / 



x (tr[(V k y] tr[(V k y)) s . (36) 

In reference j3], it was argued that in the limit I — > oo the terms with s ^ become 
negligibly small in comparison with the terms with s = 0. The latter correspond 
to the unlinked contributions. Upon resummation, these terms factorize into the 
product g(zi) ■ g(z 2 ). Therefore, the two-point correlation function R 2 (zi, z 2 ) defined in 
equation (|2*Hj) vanishes asymptotically in the limit I — > oo and, thus, implies a Poissonian 
spectrum. In detail, the argument is the following. 
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Observing that the binomial factors and the powers of Z\ and z 2 in equation do 
not depend on I or on m, and recalling that \z\\ ~ \z 2 \, we see that the only /-dependent 
factor is 



Ts = ^- 2 (tr[(V k y]tr[(V k y]) s . (37) 

Obviously To = 1 and, thus, all unlinked contributions are constant and independent 
of I. Using counting arguments, the dilute limit, and Stirling's formula, one obtains 
the asymptotic estimate T s ~ l~ sk for T s This result implies that all connected 
contributions vanish in the limit I — > 00 at least as l~ k . We note that in the dilute limit, 
l~ k vanishes as a power of the logarithm of the Hilbert space dimension N. Hence the 
two-point function R 2 (zi, z 2 ) vanishes. This line of reasoning can also be applied to the 
n-point correlation function with n > 2 showing that the spectrum is Poissonian. 

This argument has the obvious flaw that it does not display the role of e, the 
difference of the energy arguments of the two Green functions. By the same token, the 
argument does not display the difference between the connected part of the two-point 
function (where Z\ and z 2 carry infinitesimal imaginary increments of opposite signs) and 
the product of two advanced (or two retarded) Green functions (where z\ and z 2 carry 
infinitesimal imaginary increments of the same sign). In the case of canonical RMT, we 
know that in the first case {z\ and z 2 carry infinitesimal imaginary increments of opposite 
signs) the analogue of R 2 does not vanish asymptotically (in fact, the asymptotic value 
of this function determines the A3-statistic and the stiffness of the spectrum) while in 
the second case (z± and z 2 carry infinitesimal imaginary increments of the same sign) 
the analogue of R 2 vanishes asymptotically. This implies that in the case of canonical 
RMT, the analogue of R 2 has a discontinuity across the real energy axis. It is to be 
expected that such a discontinuity exists likewise for R 2 (zx, z 2 ). 

Srednicki has recently drawn attention to this singular behaviour and has 
considered the convergence properties of the sums in equation (J33j) in the light of this 
question. Equation (jHBj) is symmetric in the indices p and q. This suggests writing 
Rq(zi, z 2 ) as 

R 2 {z l)Z2 ) = t =^== T s , (38) 
where the function g s (z) is given by 



(2p-l)Uf2p + s\ 1 ^ (2 P + s)\ 
P =o z2p+s+l V s J s\z s + 1 ^ p\(2z 2 )P ■ ( ' 



In equation (|3*5|) we have used the fact that g s =o(z) = g(z), see equation (J2SJ). Using 
Borel resummation and working out the resulting integrals as in the case of the one-point 
function, one is led to the integral representation 



s+l 

(40) 



± iz u 

2 
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The upper plus (minus) sign refers to a positive (negative) infinitesimal imaginary energy 
increment, respectively. Using equation (j4*Uj) in equation (JHH|) yields 



Srednicki points out that, in order for both R 2 (zi, z£) and R<2(zf , z%) to be well 
defined, F(y) must be free of singularities on the real axis, in which case the integrals 
can be performed "without further (arbitrary) regularization" . Using the case k — 1, 
he shows that F(y) does have a singularity. This then invalidates the conclusion drawn 
in reference jlj that the spectrum is Poissonian. Srednicki concludes that the spectral 
statistics of the EGE(fc) in the dilute limit remains an unsolved problem [3E]. The 
example k — 1 considered by Srednicki is all the more remarkable because in this case, 
it is clear from independent arguments that the spectrum is Poissonian for m ^> 1. 

8. Correlated Matrix Elements: The Limiting Ensembles 

Differences between the spectral fluctuations of EGE(A;) and canonical RMT may be 
attributed to the fact that in EGE(fc), the many-body matrix elements are strongly 
correlated while in RMT, they are not. The correlations are apparent from the 
definition of the ensemble, equation (@J): Matrix elements not related by symmetry which 
involve different many-body states ipk,a, may have the same value. It is instructive to 
study whether and how these correlations of the matrix elements influence the spectral 
statistics. We do so by maintaining the graphical structure of the EGE(fc). By this we 
mean that many-body matrix elements which vanish for the EGE(fc) continue to have 
the value zero. We modify only the values of the non-vanishing matrix elements in such 
a way as to display the role of the correlations. It is in this spirit that we construct two 
limiting ensembles, one with maximum and one with minimum correlations among the 
matrix elements. Both ensembles have the same graphical structure as the embedded 
ensembles. The EGE(fc) will be seen to lie between these two limiting ensembles. For 
simplicity we consider only Fermions. We essentially follow reference 

The graphical structure of the embedded ensembles may be displayed by assigning 
to each Hilbert space vector ipj^a a vertex a, and to each non-diagonal matrix element 
(/i|Vfc|z/) which is not identically equal to zero, a link connecting the vertices fi and v. 
The diagonal matrix elements {fi\Vk\fi) are represented by loops attached to the vertex 
fi. With the exception of the loops, the resulting structure is referred to as a "regular 
graph" in the mathematical literature pH] . (For Bosons, the resulting structure is not 
a regular graph.) The number of vertices is obviously given by N, the dimension of 
Hilbert space. For Fermions, the number M of links emanating from any given vertex 



9(4) -9(4) Jo 
where the function F(y) is defined as 




(41) 



(42) 
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is the same for all vertices and given by 

M-sC)f-»)- 

This result is obtained by counting the non-diagonal matrix elements that 
connect states which differ in the occupation numbers of at most k particles. For k < m, 
we have M < N — 1 while M = N — 1 for k = m. The number P of independent links is 
given by the number of matrix elements above the main diagonal which do not vanish 
identically. Hence, 

P = ^MN. (44) 

Obviously, N, M and P do not depend on the symmetry parameter /3. 

To obtain a measure for the correlations between matrix elements, we define the 
number Kp of independent random variables. This number differs for the unitary and 
the orthogonal ensemble and is given by 




+ S, 



pi 



(45) 



The ratio Kp/P of the number Kp of independent random variables and the number 
P of different links serves as a measure of the correlations. For k <m, Kp/P is much 
smaller than one. It approaches the value (3{N + 5p\)/{N — 1) monotonically from 
below as k approaches m. This shows that for k <C m there are strong correlations 
between matrix elements on different links. The correlations disappear as we approach 
the canonical limit k = m. This is illustrated for several parameters k, m and I in 
figure |H 

We now define the two limiting ensembles by assigning a minimum and a maximum 
number of independent random variables to the graph structure of the embedded 
ensembles. For the EGE(fc), the parameter Kp/P lies between the values associated 
with these limiting ensembles. Obviously, the minimum number of random variables is 
one. The corresponding ensemble EGE min (/c) is defined in terms of the matrix elements 
V™ n \v) of the operator 

Vr n = ^£<A 7 . (46) 

The factor v is a real (complex) Gaussian random variable for /9 = 1 (fl — 2, 
respectively). Without loss of generality we may, however, put \v | 2 = 1, removing the 
distinction between the unitary and the orthogonal cases. The ensemble EGE max (fc) 
containing the maximum number of independent random variables is obtained by 
assigning, within the constraints imposed by symmetry, a different random variable 
v^v to each link of the graph. The matrix elements of EGE max (fc) are given by 

WT*\v) = <V £</#l,A>> • (47) 

For (3 = 1 (0 = 2), the matrix v^ v is real symmetric (complex Hermitean, respectively). 
Elements not connected by symmetry are uncorrelated Gaussian random variables with 
mean value zero and variance v^v^^ = S^ u >5 Ufl > + SpiS^S^i . 
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Figure 4. The ratio Ki/P, see equations (I44|l and (|45(l . on a logarithmic scale versus 
k/m. Panel (a): / = m/l= 1/5. Panel (b): / = m/l= 1/2. Panel (c): I = 100. Panel 
(d): m = 12. Taken from reference @|. 

The ensemble EGE m i n (fc) is fully integrable and has spectral fluctuations which are 
not of Wigner-Dyson type. For k — 1, one finds two different degenerate eigenvalues 
Ai = vl and A 2 = 0, with degeneracies ri\ = and n 2 = F° r k = m, 

the matrix representation of EGE m i n (m) in Hilbert space carries the entry v on every 
element. Diagonalization of this matrix is trivial and gives the eigenvalues Nv (non- 
degenerate) and zero ((N — l)-fold degenerate). Again, the ensemble is fully integrable, 
and the spectral fluctuations are not Wigner-Dyson. Using the supersymmetry method 
it can be shown that the spectral fluctuation properties of EGE max (fc) coincide with 
the predictions of RMT. From these facts, it is apparent that the limiting ensembles 
cover the extreme cases of a fully integrable system and a system with Wigner-Dyson 
spectral statistics. Figure |U then suggests that as k increases, the spectral fluctuations 
of both EGOE(/;;) and EGUE(A;) may undergo a gradual transition from Poissonian to 
Wigner-Dyson behaviour. 

9. Conclusions 

We return to the questions raised at the end of Section [T] We recall that both for k = m 
and in the dilute limit, Fermions and Bosons behave qualitatively similarly, while the 
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dense limit for Bosons is a special case. 

(i) What is the shape of the spectral density? Among the four questions formulated 
in the Introduction, this is the one to which we have a nearly complete answer. The 
moments method shows that there is a gradual transition from semicircular shape for 
k = m to Gaussian shape in the dilute limit. This is supported by the binary correlation 
approximation which yields a Gaussian shape in the dilute limit. The eigenvalue 
expansion of the matrix of second moments has added the insight that the semicircle 
prevails as long as 2k > m, and that the transition to Gaussian shape sets in at 2k = m. 
The special role of the value 2k = m is due to duality. The case of dense Bosons is 
special: It is not ergodic. 

(ii) What are the spectral fluctuation properties? The supersymmetry method 
suggests that the spectral fluctuations are of Wigner-Dyson type as long as 2k > m. It 
is possible that the range of validity of Wigner-Dyson spectral statistics extends into 
the domain 2k < m although here the supersymmetry method yields also non-universal 
contributions. In the dilute limit, the situation is not clear. The straightforward 
extension of the binary correlation approximation yields Poissonian statistics for / — > oo 
but relies on manipulations which are mathematically questionable. The eigenvector 
expansion of the matrix of second moments yields expressions which change smoothly 
with k, m, I. From these expressions one would not expect a sudden transition from 
Poissonian statistics (which surely applies for k — 1) to Wigner-Dyson statistics for 
k = 2. The group-theoretical approach shows that the part of which transforms 
under U(N m ) like a multiple of the unit matrix (and, thus, leads to Poissonian level 
statistics) is largest for k -C m. The study of the correlations between many-body 
matrix elements provides intuitive insight into the causes responsible for deviations 
from Wigner-Dyson statistics. Dense Boson systems are close to integrable and display 
no similarity to Wigner-Dyson spectral statistics. 

(iii) Are these properties universal? In analogy to canonical RMT, we ask here 
whether the spectral density and the spectral fluctuation properties hold irrespective of 
the Gaussian choice for the /c-body matrix elements, i.e., the variables v^ ai defined in 
Section |21 This is a meaningful question: We recall that in canonical RMT, ensembles 
with non-Gaussian weight factors have been studied. Since for m = k, EGE(m) is 
identical to canonical RMT, non-Gaussian ensembles of EGE(/c) type can certainly be 
defined. The analytical methods reviewed in this paper all rely heavily on the Gaussian 
distribution of the ffc ;Q . 7 's. At this point in time, no analytical approach is in sight 
to answer this question, nor are we aware of any numerical simulations addressing the 
issue. Another aspect of universality is addressed in the group-theoretical analysis of the 
embedded ensembles. It shows that the embedded ensembles for Fermions and Bosons 
are part of a much wider class of embedded random-matrix ensembles defined in terms 
of irreducible representations of some compact Lie group. 

(iv) Are the spectra ergodic? Partial answers to this question come from numerical 
simulations, from the moments method, and from the eigenvalue expansion of the 
matrix of second moments. Ergodicity can be expected only in the limit of infinite 
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matrix dimension. For k = m ergodicity holds for all observables. In the dilute 
limit, the spectral density is ergodic, but non-ergodic contributions vanish very slowly 
(with inverse powers of the logarithm of the dimension of the matrices). This causes 
difficulties in the analysis of numerically simulated spectra. For dense Bosons, even 
the low moments of the spectral density are not ergodic. It seems that this is the 
first known case of a random-matrix model which is not ergodic in the limit of infinite 
matrix dimension. We are not aware of any studies addressing ergodicity of the spectral 
fluctuations. 

This review shows that after many years of work, the determination of the spectral 
density and, especially, of the spectral fluctuation properties of EGE(/c) still poses 
difficult problems. The difficulties are due to the fact that EGE(/c) lacks the symmetry 
properties of the canonical ensembles of RMT. 
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